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Urbanisation is a fundamental phenomenon whose quantitative characterisation is still inadequate. We 
report here the empirical analysis of a unique data set regarding almost 200 years of evolution of the road 
network in a large area located north of Milan (Italy). We find that urbanisation is characterised by the 
homogenisation of cell shapes, and by the stability throughout time of high-centrality roads which 
constitute the backbone of the urban structure, confirming the importance of historical paths. We show 
quantitatively that the growth of the network is governed by two elementary processes: (i) 'densification', 
corresponding to an increase in the local density of roads around existing urban centres and (ii) 
'exploration', whereby new roads trigger the spatial evolution of the urbanisation front. The empirical 
identification of such simple elementary mechanisms suggests the existence of general, simple properties of 
urbanisation and opens new directions for its modelling and quantitative description. 

Urbanisation is a fundamental process in human history 1 , and is increasingly affecting our environment 2 
and our society 3 " 5 . The fraction of the world population living in urban areas has recently grown beyond 
50% and is expected to rapidly increase in the near future 6 . The challenges presented by such a rapidly 
urbanising world are complex and difficult to handle. They range from the increasing dependence of our society 
on fossil fuels 7 , to the emergence of socio-spatial inequalities 8 ' 9 , and to the up-rise of serious environmental 
issues 10 . Controlling urbanisation has already proven to be a difficult task, and it will become even more difficult 
in the near future, due to the large and unprecedented magnitude of urban expansion 11 , and to the intrinsic 
complexity of this phenomenon, which generates cities that are continuously changing over time 12 in a non- 
homogeneous fashion 13 . Trying to understand the elementary spatial mechanisms that govern urbanisation, 
leaving out specific historical, geographical, social and cultural factors, is nowadays more important than ever, 
especially because policy makers, professionals and researchers are actively looking for new paradigms in urban 
planning, land management and ecology 14 " 16 . The existence of acultural and non- demographic drivers of urba- 
nisation has been investigated using different approaches 17 , including classical studies in regional sciences 18 , 
theory of complex systems 19 " 21 , urban theory 22 and remote sensing 23 . A consistent amount of literature in urban 
history and morphology indicates that roads are a fundamental driver in urban evolution and, at the same time, 
one of the long-lasting constituent elements of urban forms 24 ' 25 . In the last decade, complex networks theory has 
provided significant contributions to the quantitative characterisation of urban street patterns 27,43,45 . Several 
studies have shown that road networks not only play a central role in the spatial organisation of urban areas 28,29 , 
but are important for the dynamical processes occurring on them 27,30 " 33 and for the evolution of urban systems in 
general 34,44 . However, a quantitative analysis of the historical development of urbanisation in metropolitan areas 
is still missing, and empirical evidence of the basic mechanisms governing urbanisation dynamics is still lacking. 
Our paper addresses this problem, and provides a study of the evolution over two centuries of the road network in 
a large urban area located north of Milan (Italy). 

The area under study, known as Groane, covers a surface of 125 km 2 , includes 29 urban centres within 14 
municipalities, and has essentially evolved along two main radial paths, connecting Milan to Como and Milan to 
Varese, respectively. The former path was constructed by the Romans during the II century B.C., while the latter 
was created during the XVI century. In the last two centuries, the Groane area has faced a complex process of 
conurbation, changing from a polycentric region into a completely urbanised area. This conurbation process is 
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Table 1 
Date 


List of geographical information sources used to construct the 
Source 


data set 
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Printer 


1914 


Map of Italy 


Italian Military Geographic Institute 


Raster 


1933 


Map of Italy 


Italian Military Geographic Institute 


Raster 


1955 


Aerial Photography Survey 


Italian Military Geographic Institute 


Raster 


1980 


Lombardy Regional Map 


Lombardy Region 


Raster 


1994 


Lombardy Regional Map 


Lombardy Region 


Raster 


2007 


Mosaic of Urban Municipalities Plans 


Lombardy Region 


Vectorial 



common to many large European metropolitan regions 35 . Despite 
local differences, it is possible to identify four distinct phases that 
characterize urbanization in the Groane area: i) Rural phase, (1800- 
1918): fundamentally pre-industrial, fully based on agricultural eco- 
nomy, with no major transportation infrastructures present, ii) 
Early-urban phase (1918-1945): between world wars period, still 
significantly based on agriculture economy, witnessing the first 
appearance of rail network, first small-scale sparse industrial coloni- 
sation, and limited expansion of rural settlements around the histor- 
ical centres, iii) Urban-industrial phase (1945-1990): remarkable 
sprawled residential and industrial development (especially textile 
and mechanics) along with population growth and highway con- 
struction, iv) Metropolitan post-industrial phase (1990-2012): 
decline of industrial activities, slower-paced urban sprawl, former 
polycentric organisation overwhelmed by the metropolitan con- 
tinuum caused by the merging of expanded centres, increased 
long-range mobility due to the development of high speed trains 
and large highway systems 36 . It is important to note that this area 
has never been subjected to overall large-scale planning efforts, one 
reason being that 14 different administrative bodies preside over 14 
different municipalities. By importing historical topographical and 
photogrammetrical data (for details see the Methods section) into a 
Geographical Information System (GIS) environment, we recon- 
structed the detailed road system (including minor streets) at seven 



different points in time, t = 1, 2, . . ., 7, respectively corresponding to 
years: 1833, 1914, 1933, 1955, 1980, 1994, 2007. For each time, we 
constructed the associated primal graph 27 , i.e. the graph where the 
nodes represent street junctions and the links correspond to road 
segments. The geographical information sources used to construct 
the primal graphs are reported in Table 1. In Fig. la we display the 
road networks at different times, showing how the initial small sepa- 
rate villages have grown by the addition of new nodes and links, 
eventually merging together in an homogeneous pattern of streets. 
We denote by G t = G(V t , E t ) the graph at time f, where V t and E t are 
respectively the set of nodes and the set of links at time t. The number 
of nodes at time t is then N(t) = \ V t \, while E(t) = \E t \ is the number 
of links. By definition, we have V t = V t _i UAV t and E t = E t _i U AE t > 
where AV t and AE t are, respectively, the set of new nodes and the set 
of new links added to the system in the time window ] t - 1, t] . In the 
following we will study the structure of the graph G t at different times 
U focusing in particular on the properties of the new links. As we will 
show, a quantitative analysis of the temporal evolution of the system 
will be able to reveal in particular the existence of two simple 
mechanisms that drive the evolution of the road network over time. 

Results 

Characterising the growth of the road network. The Groane area is 
basically characterised by an uninterrupted growth in the period 




Figure 1 | (a) Evolution of the road network from 1833 to 2007. For each map we show in grey all the nodes and links already existing in the previous 
snapshot of the network, and in colors the new links added in the time window under consideration, (b) Map showing the location of the Groane area in 
the metropolitan region of Milan, (c) Time evolution of the total number of nodes AT in the network and of the total population in the area (obtained from 
census data). 
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under consideration and displays remarkable modifications of the 
road system. As shown in Fig. lc, in less than two centuries, the total 
number of nodes N has increased by a factor of twenty, from the 
original 255 nodes present at t = 1 (year 1833) to more than 5000 
nodes at t = 7 (year 2007). However, the rate of growth is not 
constant over time: it is slow from 1833 to 1933, fast from 1933 to 
1980, and slow again from 1980 to 2007. These different growth rates 
are the signature of distinct phases of the urbanisation process, and 
are strongly related to the growth of the population. We indeed show 
that the number of nodes AT is a linear function of the number of 
people living in the Groane area (see Fig. 2a) or, in other words, that 
the average number of people per road intersection in the area 
remains constant over time. These results highlight the peculiarity 
of an area where the population remains almost uniformly distri- 
buted in space, and no single urban center stands out over the others. 
In order to remove the demographic component and focus solely on 
the road network evolution, we adopt the number of nodes N as the 
natural internal clock of the system, and we study the change of the 
network properties as a function of N. The number of links E grows 
almost linearly with N (Fig. 2b, top), showing that the average node 
degree (k)(t) = 2E(t)/N(t) is roughly constant over time, except for a 
slight increase from (k) ~2.57 to 2.8 when going from 1914 to 1980 
(see supplementary figure 1). In the bottom panel of Fig. 2b we report 
the total network length L tot (t)= ^2 eeE £(e) (where e denotes a 
network edge, and £{e) its length) which increases over time as 

Ltot-N? (1) 

where y ~ 0.54. Consequently, the average length of links decreases as 
iV 7-1 (see supplementary figure 1). This result is consistent with the 
evolution of almost regular two-dimensional lattices with a peaked 
link length distribution 27 ' 37 . Indeed, if we consider a uniform 
distribution of nodes with density p on a two-dimensional surface, 
the typical link length will be of order I' ~ 1 / y/p, which implies that 
the total length grows as L tot ~ E£ f ~ \/N, thus giving a value y = 1/2. 

Additional information on the structure of the road network can 
be obtained by looking at the quantity 38 

iVi+iV 3 



r N -- 



(2) 



where N k denotes the number of nodes of degree k (notice that we do 
not take into account nodes having k = 2 in the sum, since these are 
not usually considered as proper junctions). This quantity r N mea- 
sures the relative abundance of dead ends (corresponding to Ni) or 
T-shaped intersections (A/3), so that a small value of r N indicates a 
dominance of k = 4 junctions and reveals the presence of a large 



amount of grid-like patterns. Conversely, the value of r N is closer to 1 
if the network has numerous T-shaped crossings and dead ends. The 
plot of r N versus N (Fig. 2c) displays a steady decrease from — 0.87 
at year 1833 to Vm — 0.835 at year 2007, consistent with an in- 
crease of the percentage of four- ways crossings from N 4 /N~ 11% 
at year 1833 to N 4 /N~ 15.5% at year 2007. Recent studies have 
shown that the abundance of T-shaped crossings seems to be typical 
of self-organised or 'organic' urban networks, like those of Venice or 
Cairo 29 , while the grid-like layout is particular to cities whose shapes 
are the result of large-scale top-down planning efforts 38 , like 
Barcelona or New York. Consequently, one would conclude that 
what we are observing here is the evolution from an initial self- 
organised system to a rationally planned urban network. However, 
the Groane area at t = 1833 was actually a fully rural, not urban- 
organic network, and the successive evolution up to year 2007 has 
never witnessed any large-scale planning whatsoever. We therefore 
interpret the result of Fig. 2c as the signature of an evolution from a 
pre-urban' condition, with the dominance of dead-ends and 3-junc- 
tions typical of rural centres in the very early stages of growth (i.e. still 
constrained by the radial convergence of major roads), to an increas- 
ingly mature 'urban' state, in which the network expands on a rela- 
tively unconstrained land, further away from the dense radial system 
of original centres. At that later urban stage the dominance of the 4- 
junctions grid-like pattern does not result from large-scale planning, 
but from a piecemeal urbanisation, made of a multitude of small- 
scale, plot-based, scarcely coordinated developments in time, i.e. a 
substantially self- organised territorial order. 

Evolution of cells: towards homogenisation. Road networks are 
planar graphs consisting of a series of land cells surrounded by 
street segments. The statistics on the area and the shape of cells 
can be used to distinguish regular lattices from heterogeneous pat- 
terns. In particular, it has been recently observed that for Dresden 
(Germany) 28 and for a simple model of road networks 37 , the cell area 
distribution P(A) is a power law 

P(A)~A~ T (3) 

with an exponent t very close to the value 1.9. This value of z can be 
explained in terms of a lattice constructed on a set of nodes with 
density fluctuations, as shown in 27 . In Fig. 3a we report the distri- 
bution of the cell size in the Groane area at t = 2007, which is indeed a 
power law with the same exponent z = 1.9 ± 0.1. However, the 
exponent changes in time, as reported in the inset: it takes a value 
t — 1 .2 at year 1833 and converges towards t~ 1.9 as the network 
grows. Because a larger exponent indicates a higher homogeneity of 
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Figure 2 | (a) Number of nodes AT versus total population (continuous line with circles) and its linear fit (red dashed line) . (b) Total number of links E and 
total network length L toP as a function of the number of nodes N. The total network length increases as JSP 34 . (c) Value of the ratio r N between the number 
of nodes with degree k = 1 and k = 3, and the total number of nodes. In the inset we report the percentage of nodes having degree k = 4 as a function of N. 
Notice that the relative abundance of four-ways crossings increases by 5% in two centuries. 
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Figure 3 | (a) The size distribution of cell areas at t = 2007 can be fitted with a power-law p( A) ~ A~ T , with an exponent i ~ 1 .9. The inset shows the value 
of t at different times, (b) Distribution of cell shapes at different times, as quantified by the shape factor Q. The shape factor of different polygons is 
reported at the bottom axis for comparison, (c) and (d) Maps showing the cell shapes (white lines) for the network as it is before 1955 (left panel) and as it 
is after 1955 (right panel). We see on the left panel that we have predominantly triangles and rectangles, while on the right panel we can observe a 
predominance of rectangles with sides of almost the same length. 



cell areas, we are thus witnessing here a process of homogenisation of 
the size of cells. Accordingly, the relative dispersion 3 A = a(A)/fi(A) 
of cell areas (where fi(A) and a (A) are the average and the standard 
deviation of A, respectively) decreases from 0.5 at year 1833 to 0.26 at 
year 2007, indicating that the variance of the distribution becomes 
smaller as AT increases (see supplementary figure 2). 

The diversity of the cell shapes can be quantitatively characterised 
by the so-called shape factor <I> = A/(7iD 2 /4), defined as the ratio 
between the area A of a cell and the area of the circle of diameter D 
circumscribed to the cell 26 . The value of the shape factor is in general 
higher for regular convex polygons, and tends to 1 when the number 
of sides in the polygon increases. The distributions P(<I>) reported in 
Fig. 3b clearly reveal the existence of two different regimes: before 
1933 the distributions are single-peaked and well approximated by a 
single Gaussian function with an average of about 0.5 and a standard 
deviation of 0.25. Conversely, after 1955 the distributions of the 



shape factor display two peaks and can be fitted by the sum of two 
Gaussian functions. The first peak coincides roughly with the one 
observed before 1933, while the second peak, centred at 0.62, signals 
the appearance after 1955 of an important fraction of regular shapes 
such as rectangles with sides of similar lengths. In Fig. 3c-d we show 
the cell shapes at different times (before and after 1955) which visu- 
ally confirm the findings of Fig. 3. The decrease in the relative dis- 
persion 3 A and the increase in the fraction of regular shapes suggest 
that the network undergoes an evolution towards homogenisation. 
This appears to result from a combined process that exhibits two 
clear patterns: i) the fragmentation of larger cells of natural land into 
smaller ones, then heading towards urbanisation by medium-large 
manufactures or services, and ii) the mostly residential urbanisation 
of peri-urban natural land in successive rings around the historical 
urban centres. In particular, this latter pattern explains the emer- 
gence of more regularly shaped cells after the Second World War. As 
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urbanization took place increasingly further from historical main 
roads, residential blocks became less constrained by the triangular 
shape defined by those roads and, as a result, increasingly followed a 
regular rectangular shape. This regular shape is in fact the most effi- 
cient way of subdividing land for urbanization, reflecting the inner 
organization of the block into equally regular rectangular plots. This 
grid-like layout tends to be applied extensively when local constraints, 
like main roads converging into a village, do not force the development 
along different patterns (for example T-shaped junction dominated) 39 . 

Properties of new links: elementary processes of urbanization. 

Road networks grow by the addition of new streets (links) and new 
junctions (nodes). We focus here on the properties of the new links 
by looking at their length and centrality. In Fig. 4 we show the 
cumulative distribution of the length of new links according to the 
time-section in which they appeared first. The inset shows that the 
average length of new links steadily decreases over time, as expected 
from the general considerations of land fragmentation reported 
above. More precisely, we consider, at each time t, the length value 
^90%(t) such that 90% of new links at time t are shorter than € 90 %(f) 
(i.e. such that P(€ < € 90 o /o ) = 0.9). We notice that the value € 90 o /o 
decreases in the period 1833-1933 from 625 meters down to 325 
meters, while no sensible variation is observed from 1933 to 1994, 
even if the network keeps growing. In the last period, i.e. from 1994 to 
2007, we observe another decrease of € 90 %(f) from 325 to 225 meters. 
In addition, the relative dispersion of the length of the new links is 
almost constant and of order one, and the distribution does not vary 
too much after 1955. This phasing fits well with the historical 
development outlined in the introduction where during the rural 
and early-urban phase up to the second world war we observe the 
passage to an urban state that is then maintained along the core 
urbanisation age in the urban-industrial period until the 1980s, 
followed by a different post-industrial, metropolitan regime. 

The nature of the growth process can be quantitatively charac- 
terised by looking at the centrality of streets. Among the various 
centrality indices available for spatial networks, we use here the 
betweenness centrality (BC) 40-42 , which is one of the measures of 
centrality commonly adopted to quantify the importance of a node 
or a link in a graph. Given the graph G t = G( V t , E t ) at time t, the BC of 
a link e is defined as 



He)- 



(4) 



where is the number of shortest paths from node i to node j, while 
(?ij(e) is the number of such shortest paths which contain the link e. 
The quantity b(e) essentially measures the number of times the link e 
is used in the shortest paths connecting any pair of nodes in the 
network, and is thus a measure of the contribution of a link in the 
organisation of flows in the network. In order to evaluate the impact 
of a new link on the overall distribution of the betweenness centrality 
in the graph at time t, we first compute the average betweenness 
centrality of all the links of G t as 



b(G t )- 



(N(t)-l)(N(t) 



' eeE t 



(5) 



where b(e) is the betweenness centrality of the edge e in the graph G t . 
Then, for each link e* e AE t , i.e. for each newly added link in the 
time window ] t— l,f] we consider the new graph obtained by remov- 
ing the link e* from G t and we denote this graph as G t \{e*}. We 
compute again the average edge betweenness centrality, this time 
for the graph G t \{e*}. Finally, the impact db( e *) of edge e* on the 
betweenness centrality of the network at time t is defined as 



S b (e*)-- 



\b(G t ) -b(G t \{e*})] 
b(G t ) 



(6) 



The BC impact is thus the relative variation of the graph average 
betweenness due to the removal of the link e*. We have measured this 
quantity for the different times -section and we report the results in 
Fig. 5. Remarkably, the distribution of db(e) displays two well-sepa- 
rated peaks (Fig. 5d). The importance of the first peak tends to 
increase in time while the second peak decreases, until they mostly 
merge into one single peak in the last time-section (1994-2007). In 
order to understand the nature and the evolution of the two peaks, we 
focus on the geographical location of new links with different BC 
impact. In Fig. 5a-c we report respectively, the map of the network at 
year 1914, year 1994 and year 2007. We have coloured in green the 
links whose centrality impact belongs to the left peak, and in red the 
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Figure 4 | Cumulative distributions of the length of links added at different times. In the inset we report the average length (€) of new links. 
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Figure 5 | The two phases of densification (green) and exploration (red), illustrated for the networks at year 1914 (a), 1994 (b) and 2007 (c). Panel (d) 
shows the probability distribution of the BC impact d b (e) for the different time snapshots. The red peak corresponds to exploration, and the green peak to 
densification. Notice that the red peak becomes smaller and smaller with time, and completely disappears in the last snapshot, (e) Stacked chart showing 
the proportion of new links with a certain value of k min for each value of by in the network at t = 1955. The peak corresponding to exploration with a high 
BC impact is entirely due to links having k min = 1. 



links whose by corresponds to the right peak. The new links are 
divided into two classes: green links (small by, left peak) tend to 
bridge already existing streets, while red links (large Sy, right peak) 
are usually dead-ends edges branching out of existing links, generat- 
ing a new crossing and splitting the original link into two road seg- 
ments. In order to better characterise the two classes of edges we 
compute k min (e), i.e. the minimum of the degrees of both endpoints 
of a new link e, and we then consider for each value of Sy, the relative 
abundance of new links for which k min is equal to 1,2, 3, 4, respect- 
ively. We plot the results in a gray-scale stacked chart (see Fig. 5e) for 
the network at year 1955 (the plots for the other networks are shown 
in the supplementary figure 3). The analysis of k min confirms that the 
links belonging to different peaks have distinct qualitative features: 
links in the right peak have k min = 1, while links in the left peak have 
k min > 2. The distribution of BC impact thus suggests that the evolu- 
tion of the road network is essentially characterised by two distinct, 
concurrent processes: one of densification (green links, left peak, i.e. 
lower impact on centrality, k min > 2) which is responsible for the 
increase of the local density of the urban texture, and one of explora- 
tion (red links, right peak, higher impact on centrality, k min = 1) 
which corresponds to the expansion of the network towards prev- 
iously non-urbanised areas. At a closer sight, these two patterns tend 
in many cases to appear in a temporal sequence, the former acting in 
preparation for the latter, i.e. exploration being the first phase of a 



cycle of urbanisation then completed by a second phase of densifica- 
tion. Obviously, since the amount of available land decreases over 
time, at earlier time-sections (such as in the year 1914) the fraction of 
exploration is higher, while in the 80's it becomes smaller until it 
almost disappears in 2007. The only remaining peak mostly corre- 
sponds to densification with new links having k min > 2 (see supple- 
mentary figure 3). 

We finally focus on another aspect of the evolution of street cent- 
rality, namely the relation between the age of a street and its cent- 
rality. In Fig. 6a we display the links of the network in 2007 with a 
color code depending on their age, and in Fig. 6b we report their BC. 
A simple visual inspection shows that highly central links usually are 
also the oldest ones. In particular, the links constructed before 1833 
have a much higher centrality than those added at later time- 
sections. More precisely, the seven curves in panel Fig. 6c report 
the cumulative distribution of the BC computed on the network in 
2007 for the links added at the different times. This is defined as the 
probability P(b < x) that a link, appearing at a certain time- section, 
has a value of betweenness centrality b smaller than or equal to x in 
the final network in 2007. We notice that the historical structure of 
the oldest roads mostly coincide with the highly central links at year 
2007. In particular, the inset of Fig. 6c indicates that more than 90% 
of the 100 most central links in 2007 (and almost 60% of the top 
1000) were already present in 1833. This result reveals that a "back- 
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Figure 6 | Color maps indicating (a) the time of creation of each link and (b) its value of betweenness centrality (BC) at year 2007. (c) The cumulative 
distribution of BC of links added at different times. The inset reports the percentage of edges added at a certain time which are ranked in the top n 
positions according to the BC. Different curves correspond to n = 100, 500, 1000. 



bone" of highly-central routes, that have framed the Groane area in 
the rural economy of the pre-industrial period, has been driving the 
development of the area across two centuries of industrialisation, 
urbanisation and de-industrialisation up to present days, without 
any major modification. 

Urban morphologists have long observed that streets tend to per- 
sist in time much longer than other urban elements like buildings or 
land-uses 47 . Our results confirm that at another level, that of the 
topology of the street network, the persistence of roads over time 
appears to be due to their multifaceted impact on the urban form at 
different scales by deeply informing factors like accessibility to space 
and resources, land ownership and land values, each of these invol- 
ving different social actors and interests. 

Discussion 

In this paper we have studied the evolution, over almost two centur- 
ies, of the street network in a large large area close to Milan (Italy). 
Such an area is of interest for urban studies because it displays an 
important process of urbanisation. Urbanisation is reflected in the 
growth of the road network and displays different speeds at different 
times, with a fast increase occurring between 1933 and 1994. Our 
results reveal a quantitative signature of urbanisation on the evolu- 
tion of the shape and size of land cells, which become more homo- 
geneously distributed and square- shaped. Simultaneously, we 
observe a general trend towards a larger number of 4-way junctions, 
as opposed to an earlier structure of dead ends and 3 -way junctions. 
These structural transformations appear to be the result of the inter- 
play between two concurrent dynamics, namely densification and 
exploration. While exploration is typical of the earliest historical 
periods of urbanisation, densification predominates in the latest. 

We were also able to quantitatively characterise the stability of the 
structure of most central streets over time. For instance, the most 
central streets in the network in 2007 largely coincide with the oldest 
ones. Central roads appear to therefore constitute a robust spatial 
backbone which remains stable over time, and characterises the 
evolution of the road system as a continuous expansion and rein- 
forcement of pre-existing structures rather than as a sharp switch 
towards radically new configurations. The kind of evolution that we 
are witnessing in the Groane area is therefore certainly massive in 



overall quantitative terms, reflecting unparalleled changes in the pace 
and nature of economic, technological and social order during the 
transition from a mainly agricultural to a post-industrial age; how- 
ever, it is still an organic form of evolution that builds at each step on 
the achievements of the previous, confirming in the long term a 
resilient structure that tends to persist in time all over the urbanisa- 
tion process. These findings quantitatively support the hypothesis 
that spatial systems undergoing fine-grained forms of evolution tend 
to exhibit simple local dynamics of change, continuously expanding 
upon pre-existing structures, a morphogenetic behaviour that they 
appear to share with many living organisms and other natural 
systems 46 . 

The results of this research, including the existence of the two 
main dynamics of exploration and densification, together with the 
persistence over time of a structural backbone made of highly central 
streets, cannot be extended to the generality of urbanisation pro- 
cesses unless supported by further investigation of cases in different 
geographical and economic positions. An important point in our 
view is that this study provides quantitative results which can be 
confirmed (or falsified), a critical aspect in Science. 

Methods 

Temporal Network Data. The area under study covers 125 km 2 and includes 29 
urban centers within 14 municipalities that have developed along two main radial 
paths, connecting Milan to Como and Milan to Varese. We sampled seven street 
networks primal graphs at different times, where the street junctions are represented 
as nodes and the roads (or streets) are the links of the networks. We used a mixed 
procedure of ArchMap tool extensions and python scripts operating over a geo- 
database. First, the street network has been drawn in vectorial format on the basis of a 
collection of historical areal images and historical maps imported in the ArcGISlO 
platform. Sources used for each temporal steps are reported in Table 1. The images 
have been imported as a .geotif extension, and all the layers has been projected using 
the coordinate system Monte Mario (Rome), Italy zone and Transverse Mercator 
projection. The first geo- referencing (alignment of different images) has been 
performed on the basis of historical buildings and landmarks such as paths and roads. 
The redraw operation was done using an ad-hoc python script for creating the node 
layers representing the street junctions. The second geo -referencing operation has 
been performed over the street junctions using the spatial analyst tool with the spatial 
join utilities provided in ArchMap 10. Subsequently, for all time layers, using an ad- 
hoc python script we produced the weighted adjacency lists that has been used for the 
network analysis, the weight correspond to the real length of the street. We denote by 
G t = G( V t , E t ) the obtained primal graph at time t, where V t and E t are respectively the 
set of nodes and links at time t. The number of nodes at time t is then N(t) = | V t \ and 
the number of links is E(t) = \E t \. Using common definitions, we thus have 
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V t = Vt-iKjAVt and E t = E t -i{JAE t , where AV t and AE t are respectively the new 
street junctions and the new streets added in time ]t - 1, t] to the network existing at 
time t-l. 



1. Mumford, L. The City in History (Harcourt Brace, New York, 1961). 

2. McDonald, R. L, Green, P., Balk, D., Fekete, B. M., Revenga, C, Todd, M. & 
Montgomery, M. Urban growth, climate change, and freshwater availability. Proc. 
Natl Acad. Sci. (USA), 108, 6312-6317 (2011). 

3. Dye, C. Health and Urban Living. Science 319, 766-769 (2008). 

4. Glaeser, E. Cities, Productivity, and Quality of Life. Science 333, 592-594 (201 1). 

5. Barnes, I., Duda, A., Pybus, O. G. & Thomas, M. G. Ancient urbanization predicts 
genetic resistance to tuberculosis. Evolution 65, 842-848 (2011). 

6. UN Population division, http://www.unpopulation.org. 

7. Newman, P. & Kenworthy, J. Sustainability and Cities: Overcoming Automobile 
Dependence (Island, Washington, DC, 1999). 

8. UN-Habitat. The challenge of slum, Global Report on Human Settlements 
(Eacrthscan, London and Sterling, VA, 2003). 

9. UN-Habitat. State of world's cities 2010/11, bridging the urban divide (Earthscan, 
London and Sterling, VA, 2008). 

10. Hardoy, J. E., Mitlin, D. & Satterthwaite, D. Environmental Problems in an 
Urbanizing World (Earthscan, London, 2001). 

11. Angel,S., Sheppard, S. C. & Civco, D. L. The Dynamics of Global Urban Expansion 
(The World Bank, Washington, DC, 2005). 

12. Hall, P. Modelling the Post-industrial city. Future 294-5, 311-322 (1997). 

13. Seto, K. C, Sanchez-Rodriguez, R. & Fragkias, M. The New Geography of 
Contemporary Urbanization and the Environment Annu. Rev. Environ. Resour. 
35, 4.14.28 (2010). 

14. UN-Habitat. Planning sustainable cities. (United Nations Human Settlements 
Programme, Nairobi, Kenya, 2010). 

15. Turner, B. L., II, Lambin, E. F. & Reenberg, A. The emergence of land change 
science for global environmental change and sustainability. Proc. Natl Acad. Sci. 
(USA), 104, 20666-20671 (2007). 

16. Ramalho, C. E. & Hobbs, R. J. Time for a change: dynamic urban ecology. Trends 
in Ecology & Evolution 27, 179-188 (2012). 

17. Batty, M. The Size, Scale, and Shape of Cities. Science 319, 769-771 (2008). 

18. Wilson, A. G. Entropy in Urban and Regional Modelling (Pion Press, London, 
1970). 

19. Batty, M. & Longley, P. A. Fractal Cities: A Geometry of Form and Function 
(Academic Press, San Diego, CA, 1994). 

20. Batty, M. Cities and complexity (The MIT Press, Cambridge, MA, 2005). 

21. Makse, H. A., Havlin, S. & Stanley, H. E. Modelling urban growth patterns. Nature 
377, 608-612 (1995). 

22. Salingaros, N. Principles of Urban Structure (Techne Press, Amsterdam, Holland, 
2005). 

23. Taubenbck, H. Esch, T., Felbier, A., Wiesner, M., Roth, A. & Dech, S. Monitoring 
urbanization in mega cities from space. Remote Sensing of Environment 117, 
162176 (2012). 

24. Marshall, S. Streets and Patterns (Spon Press, Abingdon, Oxon UK, 2006). 

25. Southworth, M. & Ben-Joseph, E. Streets and the Shaping of Towns and Cities 
(Island Press, Washington DC. USA, 2003). 

26. Haggett, P. & Chorley, R. J. Network analysis in geography. (Edward Arnold, 
London, 1969). 

27. Barthelemy, M. Spatial Networks. Physics Reports, 499, 1-101 (2011). 

28. Lammer, S., Gehlsen, B. & Helbing, D. Scaling laws in the spatial structure of 
urban road networks. Physica A, 363, 89-95 (2006). 

29. Cardillo, A., Scellato, S., Latora, V. & Porta, S. Structural properties of planar 
graphs of urban street patterns. Phys. Rev. E 73, 066107 (2006). 

30. Bettencourt, L. M. A., Lobo, J., Helbing, D., Kuehnert, C. & West, G. B. Growth, 
innovation, scaling, and the pace of life in cities. Proc. Natl Acad. Sci. (USA), 104, 
7301-7306 (2007). 

31. Bettencourt, L. M. A. & West, G. B. A unified theory of urban living. Nature, 467, 
912-913 (2010). 



32. Balcan, D., Colizza, V., Goncalves, B., Hu, H., Ramasco, J. & Vespignani, A. 
Multiscale mobility networks and the spatial spreading of infectious diseases. 
Proc. Natl Acad. Sci. (USA), 106, 21484-21489 (2009). 

33. Porta, S., Latora, V., Wang, F., Rueda, S., Strano, E., Scellato, S., Cardillo, A., Belli, 
E., Cardenas, F., Cormenzana, B. & Latora, L. Street Centrality and Location of 
Economic Activities in Barcelona. Urban Studies Nov. 10 2011, 1-18 (2011). 

34. Hanson, S. & Giuliano, G. (Eds). The geography of urban transportation. (Guilford 
Publications, New York, NY, 2004). 

35. European Commission, Joint Research Centre. Urban sprawl in Europe, The 
ignored challenge. EE A Report 10, (2006). 

36. Turri, E. La Megalopoli Padana (Marsilio, Venezia, IT, 2000). 

37. Barthelemy, M. & Flammini, A. Modeling urban street patterns. Phys. Rev. Lett. 
100, 138702 (2008). 

38. Courtat, T., Gloaguen, C. & Douady, S. Mathematics and morphogenesis of cities: 
A geometrical approach. Phys. Rev. E 83, 036106 (2011). 

39. Ben- Joseph, E. The Code of the City: Standards and the Hidden Language of Place 
Making (MIT Press, Cambridge, MA, 2005). 

40. Freeman, L. C. A set of measures of centrality based on betweenness. Sociometry 
40, 35-41 (1977). 

41. Crucitti, P., Latora, V. & Porta, S. Centrality measures in spatial networks of urban 
streets. Phy. Rev. E 73, 036125 (2006). 

42. Porta, S., Crucitti, P. & Latora, V. The network analysis of urban streets: a primal 
approach. Environment and Planning B: Planning and Design 33, 705-725 (2006). 

43. Latora, V. & Marchiori, M. A measure of centrality based on network efficiency. 
New Journal of Physics 9, 188 (2007). 

44. Xie, F. & Levinson, D. Topological evolution of surface transportation networks. 
Computers, Environment and Urban Systems 33, 211-223 (2009). 

45. Xie, F. & Levinson, D. Measuring the structure of road networks. Geographical 
Analysis 39, 336-356 (2007). 

46. Alexander, C. The Nature of Order: an essay on the art of building and the nature of 
the universe. Book 2, The Process of Creating Life (Oxford University Press, 
Oxford, UK, 2003). 

47. Vernesz-Mouton, A. Built for change: neighborhood architecture in San Francisco 
(MIT Press, Cambridge, MA, 1989). 



Acknowledgements 

MB thanks R. Morris for useful comments and the Morphocity group (P. Bonnin, P. Bordin, 
C.-N. Douady, S. Douady, J.-P. Frey, P. Vincent) for stimulating discussions. ES and SP are 
grateful to the Laboratory of Territorial Planning at DIAP, Polytechnic of Milan (S. Serini, 
A. Benedetti L. Terlizzi, G. Graj and Prof. P. L. Paolillo) for providing cartographic material 
and vectorial street layers. They also thank A. Venerandi for his help in the first phase of this 
study. 

Author contributions 

All the authors have equally contributed to the design of the experiment, to the analysis and 
interpretation of the results and to the preparation of the manuscript. 

Additional information 

Supplementary Information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 
License: This work is licensed under a Creative Commons 

Attribution-NonCommercial-NoDerivative Works 3.0 Unported License. To view a copy 
of this license, visit http://creativecommons.Org/licenses/by-nc-nd/3.0/ 

How to cite this article: Strano, E., Nicosia, V., Latora, V., Porta, S. & Barthelemy, M. 
Elementary processes governing the evolution of road networks. Sci. Rep. 2, 296; 
DOI:10.1038/srep00296 (2012). 



SCIENTIFICREPORTS | 2 : 296 | DOI: 1 0.1 038/srep00296 



8 



